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The characterization of the set of quantum correlations in Bell scenarios is a problem of paramount 
importance for both the foundations of quantum mechanics and quantum information processing in 
the device-independent scenario. However, a clear-cut (physical or mathematical) characterization 
of this set remains elusive and many of its properties are still unknown. We provide here simple and 
general analytical conditions that are necessary for an arbitrary bipartite behaviour to be quantum. 
Although the conditions are not sufficient, we illustrate the strength and non-triviality of these 
conditions with a few examples. Moreover, we provide several applications of this result: we prove a 
quantitative separation of the quantum set from extremal nonlocal no-signaling behaviours in several 
general scenarios, we provide a relation to obtain Tsirelson bounds for arbitrary Bell inequalities 
and a construction of Bell expressions whose maximal quantum value is attained by a maximally 
entangled state of any given dimension. 
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I. INTRODUCTION 

Bell’s theorem shows that the correlations exhibited by 
quantum mechanical systems go beyond what is achiev¬ 
able by any local realistic theory (i. e. any local hidden 
variable model) [1]. This result has deep implications and 
still spurs several subfields of research [2], On the one 
hand, from a foundational perspective, it has been ques¬ 
tioned what physical principle in nature could then be 
responsible for giving rise to precisely the particular set 
of correlations allowed by quantum mechanics. Since it 
has been shown that this set is strictly smaller than what 
could be achieved by just imposing the no-signaling prin¬ 
ciple (i. e. the impossibility of instantaneous propagation 
of information) [3, 4], several works have considered more 
restrictive physically-motivated axioms [5-7]. Although 
these approaches rule out several subsets of no-signaling 
correlations, no definite answer has yet been found as 
there still exist supra-quantum correlations compatible 
with these principles [8]. On the other hand, from a 
more practical point of view and in the context of quan¬ 
tum information theory, it has been realized that quan¬ 
tum nonlocality can be regarded as a resource for device- 
independent quantum information processing (DIQIP). 
The tasks that can be carried out in this way include key 
distribution for cryptography [9], randomness generation 
[10] and dimensionality certification (see e. g. [11, 12]). 

In order to understand which physical principles con¬ 
strain the set of quantum correlations and to elucidate 
what are the ultimate limitations behind DIQIP pro¬ 
tocols, a fundamental question arises: can an efficient 
mathematical description of this set be found? It turns 
out that this is a highly nontrivial question. So far, the 
only systematic and general form to bound this set is 
given by the NPA hierarchy [13]. This provides a hi¬ 
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erarchy of semi-definite programs that approximate the 
quantum set from the outside. Although this constitutes 
an extremely powerful tool, it has been mainly exploited 
numerically. Thus, many properties of the quantum set 
remain unknown, since verifying general statements can¬ 
not be usually approached by numerical means. It would 
be therefore desirable to have simple analytical condi¬ 
tions constraining the set of quantum correlations. These 
could be moreover used to exclude some general subsets 
of no-signaling correlations from the quantum set based 
on their analytical properties or to bound the maximal 
efficiency DIQIP tasks with a given structure can at¬ 
tain. The aim of this article is to fill this gap: I pro¬ 
vide simple general analytical conditions any quantum 
behaviour should satisfy which rely on standard matrix 
analysis. Although these conditions emerge from the first 
step of the NPA hierarchy and are, therefore, only nec¬ 
essary, I will consider examples showing their strength 
and non-triviality. Furthermore, I will provide several 
different applications of this result: a proof of the separa¬ 
tion between the quantum set and extremal nonlocal no¬ 
signaling correlations (a question recently raised in [14]), 
a systematic way to obtain quantum bounds on arbitrary 
Bell inequalities (generalizing the recent result [15]) and 
the possibility to do robust self-testing of bipartite max¬ 
imally entangled states by building Bell inequalities that 
are maximally violated by them. 


II. FRAMEWORK AND MAIN RESULT 

We will consider the standard bipartite ( mAmsdAdB ) 
Bell scenario [2] in which two parties A and B that could 
have interacted in the past remain now uncommunicated. 
Party A ( B) can freely choose questions from a finite al¬ 
phabet X = {1,2,..., tua} (J- 7 = {1,2,..., tub})- Given 
any of these inputs x 6 X and y € y, each party can 
obtain an outcome a £ A = { 1 , 2 ,..., cU} and b £ B = 
{ 1 , 2 ,..., cIb}, where the output alphabets are also finite. 
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The central object here, referred to as behaviour, is the 
joint conditional probability distribution of obtaining the 
outputs (a, b) given the choice of inputs (x,y), P(ab\xy) 
(for which we will use the shorthand P). This list of 
dAdBxnAmB numbers must fulfill P(ab\xy) > 0 Va, b, x, y 
and b P ( ab \ x y) = 1 Va;, j/. Moreover, since communi¬ 
cation among the parties is not possible during the choice 
of input and recording of the output, the marginal of 
each party must be independent of the other’s action, 
P(a\x) = J2 b P(ab\xy) = Va,x,y ± y' and 

P(b\y) = 'Ea P ( ab \ x y) = J2a p ( ab \ x 'y ) V6,y ,x ± x'. 

The set of all behaviours satisfying these no-signalling 
conditions will be denoted by A (S. Particular elements 
of this set are the different local deterministic 

behaviours (LDBs) Di(ab\xy) = S a ,f i {x) b b,g i (y) i n which 
for every party a unique output occurs with probability 
1 for every choice of input. The convex hull of these be¬ 
haviours gives rise to the local set £ [16]. On the other 
hand, the quantum set Q is given by all behaviours that 
can be obtained by performing measurements on bipar¬ 
tite quantum states of unrestricted dimension pab, i- e. 
P(ab\xy) = tT(pABEa^Fb) f° r some projectors {E%, F b } 
such that Y a p a and Y2b p l! equal the identity in each 
party’s Hilbert space Vx, y. The crucial observation men¬ 
tioned in the introduction is that £ C Q C A fS. Al¬ 
though Q is a convex set, it is in general very hard to 
decide from the definition whether a given behaviour is 
in Q or not. On the contrary, £ and J\fS are both con¬ 
vex polytopes with vertices given by the LDBs in the first 
case, to which we have to add some non-local vertices in 
the second case. Following the standard notation, we will 
refer to these extremal non-local behaviours as Popescu- 
Rohrlich (PR) boxes. 

In order to present our results, we will use some further 
notation. We will arrange every P £ AfS to form the 
rriAdA x rnsds = tia xn^ real matrix 

p = p (ab\xy)\xa)(yb\, (1) 

abxy 


where in the standard notation of quantum mechanics 
| xa) = |x)®|a) and {|x}} denotes the computational basis 
of R mA and similarly for the other alphabet elements. 
Thus, P can be partitioned as a block matrix 


P = 


( p i 


V Pm A l 


Pl 


lm,B 


Prr 




X ns 


with blocks 


Pr 


( P(ll|xy) 
\ P(d A l\xy) 


P{ld B \xy) 


P(d A d B \xy) 


g gd^XtiB 


(column) is equal. We will furthermore consider different 
matrix norms. Following the Schatten p-norm notation, 
|| • ||i will be the trace norm (i. e. the sum of all singular 
values) while || • ||oo the spectral norm (i. e. the maximal 
singular value). We are now in the position to state our 
main result. 

Theorem 1. In every () scenario, if P £ Q 
then ||P||i < y/iriAmB■ 

We will actually prove the following stronger result: 
for every P £ Q and G £ R nAXraB it must hold that 

(P, G) = tr (PG T ) = p (ab\xy)G{ab\xy) 

abxy 

< \\G\\ 00 y/mA'm,B- ( 2 ) 

The relevance of this result will be discussed later on. 
For the moment, let us simply point out that Theorem 
1 follows from it by noticing that ||P||i = maxo tr(PO) 
where the maximization is over all isometries [17]. 

Proof. In order to prove that inequality (2) is true, we 
will show that the inequality holds VP £ Q 1 , i. e. the 
first step of the NPA hierarchy (notice that Q C Q 1 ). 
It is worth mentioning that Theorem 1 can be proved 
without invoking this [18]. However, we have chosen to 
present this proof because it furthermore clarifies the 
relative strength of the condition. Behaviours in Q 1 
must fulfill [13] that a (ha + tib) x {jia + tib) real pos¬ 
itive semidefinite matrix T = ^ p T ^ ^ exists with 

Qa = Pa(*) and Rjj = Pb (j ), where Pa ( Pb ) is a vector 
of K mAdA (M. mBdB ) with entries given by P(a\x) ( P{b\y )) 
with xa (by) in lexicographical order [42]. Thus, defining 

, for any given G the maximum value of 

(P, G) attainable in Q cannot be larger than 

maxtr(rW)/2, subject to 

r > o, 

t T(Dfr)=P A (i) (i = l,...,n A ), 
tr(Pf r ) = P B (j) (j = l,...,n B ), 

where Df = Ei ® 0„ B (Df = 0 nA ® Ef) and Ei (Ej) 
a ha x tia (n B x n B ) matrix whose only nonzero entry 
is the ii ( jj ) with value 1. This is the primal form of a 
semi-definite program (SDP) [19] with cost function p(T) 
and we will denote its solution by p(T*). The dual form 
of this SDP corresponds to 

^ , subject to diag(x) — W/2 > 0, 


t f Pa 
mm x 0 

X \ 


w = 


0 G 
G T 0 


Normalization imposes that the entries in each block add 
up to one while no-signaling that the sum of the elements 
in the same row (column) for blocks in the same row 


where x is a ua + n B real vector yielding the value d(x). 
By duality, for any feasible x (i. e. satisfying the con¬ 
straint above), it must hold that p(r*) < d(x). Thus, to 
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finish the proof it suffices to construct a feasible x yield¬ 
ing the value d(x ) = | (GHoo This is the case for 

x = \\GWoo/2(^/m B /m A l nA , \Jm A /m B l nB ) where 1„ is 
the n-dimensional vector with all entries equal to one. 
To see that it is feasible amounts to checking that 

\J m B /m A 11G11 oo In A -G 

-G T \/'m A /mB\\G\\ 00 t. nB 

This is indeed true because, since the upper left corner 
in strictly positive definite, by Schur’s complement con¬ 
dition [20] this is equivalent to 

_/ G T G \ 

yJm A /m B M|G||ooln B — || G || 'J — 0, 

which is obviously true given that the maximal eigenvalue 
of G T G is precisely || G\ |□ 

One of the appealing properties of Theorem 1 and in¬ 
equality (2) is that they have a very compact form. How¬ 
ever,the reasoning used in this proof can be applied to 
obtain other stronger but more complicated conditions 
[43]. Let us denote by M the matrix constructed using 
the same prescription as P in Eq. (1) but with entries 
given by 



M(ab\xy) = P(ab\xy) — P(a\x)P(b\y). (3) 

We will now prove that for any choice of matrix G £ 


(M,G) <||G||ooV m A^s - ^ x (^TP(a\xY 

2 V TJlA 

Y ax 

(4) 

' by 

should hold VP £ Q and, hence 

Theorem 2. In every (m A m B d A d B ) scenario, if P £ Q 
then 


||M||i < yJm A m B 


y- P(a\x) 2 

2m A 


V P(b\y) 2 
J 2m b 


Proof. This proof goes along similar lines as the previous 
one, so we will only sketch the details. The set Q 1 is 
actually equivalent to the positive semidefiniteness of T — 

( Ps ) ( P a p b ) [6], that is, 


( Q M 
\M T R 


> 0 , 


than the primal SDP 

maxtr(fW)/2, subject to 
f 

f > 0, 

tr(Dff) = P A (i) - P A (i) 2 = n A ), 

tr{Df f) = P B (j) - P B (j) 2 (j = 1,..., n B ), 
with dual 

mina; T ( D J , subject to diag(x) — W/2 > 0, 
x \4s r B J 

where the vectors P AB have entries P A (i ) 2 and P B {j) 2 . 
The same choice of x as in the previous proof does the 
job. □ 

Let us finish this section analyzing the strength of The¬ 
orems 1 and 2 by considering a few examples. As men¬ 
tioned before, our conditions are deduced from the defini¬ 
tion of the first step of the NPA hierarchy, Q 1 . Obviously 
then, the best we can expect from them is to separate 
this set from its complement in J\fS. Therefore, since 
Q C Q 1 , it is clear that there exist supra-quantum be¬ 
haviours satisfying the conditions of Theorems 1 and 2 (i. 
e. they are necessary for a behaviour to be in Q but not 
sufficient). We have performed numerical explorations 
in the (2222) scenario that show that the gap between 
Q 1 and the behaviours satisfying our conditions is small. 
Theorem 2 only provides a slightly better approximation 
of Q 1 than Theorem 1. A more detailed example can be 
found in Figure 1. To my knowledge, the only previous 
instance of analytical means to constrain Q is given by 
the results of [21] (which emerge from Q 1 as well). How¬ 
ever, the application of these techniques is not completely 
straightforward as they rely on some choice of data pro¬ 
cessing. Still, the inequality emerging from this approach 
reproduces Q 1 in the (2222) scenario. Notwithstanding, 
Figure 1 shows that our conditions give a much tighter 
restriction already in the (3322) example considered in 
[21]. This is also apparent in the (2233) case, where the 
results presented in [21] fail to completely reproduce Q 1 
for the isotropic behaviours obtained by mixtures of a 
fully random box and the PR box Pp B ( 2, 3) (see Eq. (7) 
below). On the contrary, Theorems 1 and 2 are tight in 
this case. Thus, it is interesting to point out that while 
inequalities (2) and (4) depend on a proper choice of G 
for each case, these examples indicate that Theorems 1 
and 2, which are straightforward to apply in general, still 
provide remarkably strong conditions. To show further 
the usefulness of these results, in the remaining sections 
we provide several applications of them. 


III. NONTRIVIALITY OF THE CONDITIONS 
AND NONQUANTUMNESS OF EXTREMAL 
NO-SIGNALING BEHAVIOURS 


where now Qa = P A (i) — P A {i) 2 and Rjj = P B {j ) — The previous examples already give a good idea of the 
P B (j) 2 . Thus, the value of (M,G) in Q cannot be larger strength of the conditions derived here. To analyze in 
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FIG. 1: A slice of the (3322) polytope spanned by behaviours 
P = qPi + pP 2 + (1 — p — q)P n , where Pi and P 2 are two 
PR boxes and P n is the fully mixed behaviour. Horizontal 
lines achieve the same value for the I3322 Bell inequality (see 
[21] for details). Crosses: The boundary of £ and Q (equal 
in this slice, 13322 = 0). Dashed line: The boundary of Q 1 
(I3322 = 0.2). Solid line: The conditions of Theorems 1 and 2 
(they are equivalent in this case, 0.2 < I3322 < 0.24). Dashed- 
dotted line: The inequality of [21] with the sign binning data 
processing (13322 — 0.41). A more clever data processing for 
this case only manages to lower this condition to I3322 — 0.3 
[ 21 ]. 


full generality their non-triviality, notice that, due to the 
convexity of the trace norm, ||P||i must attain its maxi¬ 
mum value in J\fS at the vertices of the polytope. Hence, 
the ideal situation would be that the LDBs achieve the 
maximal possible value (^/m^mp) and that all PR boxes 
violate this bound. Actually, it has been only recently 
shown in [14] (see also [7]) that all PR boxes (includ¬ 
ing the multipartite case) are not in Q. It was left as 
an open question there whether there exists a separation 
between them and Q. The interest of this question relies 
on the fact that PR boxes are the most advantageous re¬ 
sources in DIQIP applications such as cryptography [14]. 
Thus, the closer Q can be to these behaviours the more 
efficient these applications can be. In the following we 
show that all LDBs attain the , Juiaitib bound and that 
all PR boxes in (22d^dg) and (mm 22 ) scenarios violate 
the bound. Thus, besides showing that the bound is not 
improvable and in general not trivial, we further provide 
a simple proof in these cases of the result of [14]. More¬ 
over, we show that there actually exists a quantitative 
separation between these PR boxes and Q, answering in 
these scenarios the question raised therein. 


Proposition 3. In every ( mAmsdAds) scenario, if P £ 
C then ||P||i < y/mAfUB, with equality for LDBs. 


Proof. Every LDB is of the form Di(ab\xy) = 
df(a\x)df(b\y) where the lists df(a\x) {df(b\y)) have 


itia (tub) entries equal to 1 and 0 otherwise. Hence, 



J2 d i(. b \y)( b y\ 


1 


df(a\x)\ax) 



ax 

2 

by 


= ^ /mAmB■ (5) 


Notice that the fact that ||P||i < y/mAiriB holds VP £ C 
follows then by the convexity of the trace norm without 
the need of invoking Theorem 1. □ 


Evidently, that the inequality is fulfilled in C was al¬ 
ready obvious from Theorem 1. The important observa¬ 
tion here is that all LDBs attain the bound, hence show¬ 
ing that it cannot be improved. 

We analyze now the values ||P||i might take for PR 
boxes. Let us consider first the {f22dAds) scenario. Tak¬ 
ing Ad to be the circulant d x d matrix 


/ 0 0 ... 1 \ 


Ad = 


■■ -.o 

\ 0 ...0 1 0 / 

the PR boxes in this case are given by [44] 


( 6 ) 


Ppr^i d) 


1(1:, 

d \ 1 d Ad J ’ 


2 < d < min(d J 4 , dg), (7) 


up to relabelings of the inputs and the outputs [22]. Since 
these transformations amount to certain permutations of 
the rows or columns of P that leave the trace norm in¬ 
variant, it suffices to compute it for the matrix given in 
Eq. (7). Using the pinching inequality [23, 24], we obtain 
that 


||Pp fl (2, d)||i > - 


Id 0d 

0d A d 


= 2 . 


The conditions under which equality is attained in the 
pinching inequality are given in Theorem 8.7 of [24] and 
it is easily checked that they are not met in this case. 
Hence, we obtain that ||Ppg(2, d)| |i > 2, which amounts 
to the non-quantumness of these PR boxes by Theorem 1. 
Moreover, by a more refined use of the pinching inequal¬ 
ity, we obtain the following stronger result, which shows 
the existence of a finite gap between these PR boxes and 

Q. 


Theorem 4. In every (22d^dg) scenario, it holds VP £ 
Q that Vd \\P PR ( 2 ,d)-P \\ 1 > ||P Pfl (2,d)|| 1 - ||P||i > 
VE-2. 
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Proof. We need to show that ||Ppfl(2, d)||i > y/5 Vd. 
By permutation matrices, that leave the trace norm in¬ 
variant, we can map P = Yhabxy P{ a ^\ x y)\ xa ){v^\ to 
P' = J2 abxy P(ab\xy)\ax){by\, which has now blocks 
given by 


/ P(a6|ll) 


P' - 
ab 


P(ab\lrriB ) 


V P(ab\m A l) ■■■ P{ab\mAiTiB) 

We therefore have that 

||P P p(2,d)|| 1 = ||P^(2,d)|| 1 

( Z 0 2 ... Y\ 


Y 

0 2 

. ■■■ o 2 

V 0 2 ... 0 2 Y Z ) 


XrriB 


where 


Z = 


1 1 
1 0 


Y = 


0 0 
0 1 


( 8 ) 


Hence, using now the pinching inequality we obtain the 
desired result | |Ppp(2, d)||i > (d\\Z\\f)/d = y/§. □ 


IV. TSIRELSON BOUNDS AND BIPARTITE 
MAXIMALLY ENTANGLED STATES 

One can see that the left-hand-side of inequality (2) de¬ 
fines an arbitrary Bell expression, i. e. any linear combi¬ 
nation of the elements P(ab\xy). Since C, Q and ATS are 
compact convex sets, there always exist such expressions 
separating them, i. e. (P, G) < Gc,Gq,Gj^s depend¬ 
ing on whether P G C, Q,MS with Gc < Gq < Gj^s- 
The most characteristic one is the CHSH inequality (see 
below) in the (2222) scenario for which Gc = 2 [26], 
Gq = 2y/2 [27] and Gjys = 4 [4]. While to determine 
the optimal value of Gc and Gjyfs it suffices to check 
over the corresponding vertices, to determine the opti¬ 
mal value of Gq, known as Tsirelson bounds, is a less 
straightforward task [45]. However, this is very relevant 
to identify optimal DIQIP performances in the context 
of quantum games [2, 28]. Thus, a remarkable feature 
of inequality (2) is that it provides a systematic way to 
construct quantum upper bounds to arbitrary Bell in¬ 
equalities. Actually, our result resembles that of [15, 29], 
but the latter only holds for the particular class of Bell 
inequalities based on correlators. Later on, we will dis¬ 
cuss this relation further. To give a hint of the usefulness 
of inequality (2), we will show now that it allows to ob¬ 
tain Tsirelson’s bound for the CHSH inequality. This can 
be expressed by a matrix Gchsh with blocks 


It might be interesting to note that ||-Ppp(2,rf)|< 
2y/2. To see this, notice that for the Frobenius norm 
||Ppp(2,d)|| 2 = \/tr(Ppp(2,d)Ppp(2, d) T ) = 2/Vd and 
that ||X||i < v / n||X|| 2 for any n x n matrix X. Nu¬ 
merics suggest that the above estimates can be im¬ 
proved to 1 + y/~2 = ||Ppp(2, 2)11 1 < ||Ppp(2,d)||i < 
lim^oo ||Ppp(2,d)||i ~ 2.55, that would change the 
bound in Theorem 4 to y[2 — 1. 

Let us move now to the (mm22) scenario. The corre¬ 
sponding PR boxes have all been determined in [25] (see 
Table II therein). We denote an arbitrary one of them by 
Ppp(?n, 2). One sees that (up to relabelings) these matri¬ 
ces always have the following structure: they have a 4 x 4 
block in the diagonal given by Ppp(2, 2) followed by m—2 
2x2 blocks in the diagonal, which are either l 2 /2, A 2 /2 
or diag( 1,0). Since the latter blocks all have unit trace 
norm, it follows again by the pinching inequality that 
\\Ppii(rn, 2) ||i > ||Ppk(2, 2)||i+to — 2 = m + \/2 — 1. 
This shows again by virtue of Theorem 1 that there is 
a quantitative separation between the PR boxes in these 
scenarios and Q: 

Theorem 5 . In every ( mm22 ) scenario, it holds VP € Q 
that \\Ppn(m, 2) — P||i > y/2 — 1. 

It might be interesting to mention as well that in the 
(mm22) scenario for the fully nondeterministic boxes it 
holds that ||Ppp(m, 2)||i < my/rn. For these behaviours 
there are 2 to 2 non-vanishing entries with value 1/2 [25]. 
Hence, similarly as before, ||Ppfl(m, 2)|| 2 = m/y/ 2, ob¬ 
taining the above estimate. 


Gil — Gl2 — G 2 1 — —G 2 2 — 



It turns out that 11 Gchsh 11oo = 2 and, hence, we obtain 
the trivial bound 4. Nevertheless, given that behaviours 
in AfS must fulfill several different constraints, equivalent 
Bell inequalities can be expressed up to rescaling and 
addition of an offset. Thus, if we take (/ is a matrix with 
all entries equal to 1) 


G' f 


CHSH 


h o 2 

0 2 I 2 


+ m GcRSH ' 


it holds that (P, G' CHSH ) = 1 + (P, Gchsh)/ (2\/2) VP e 
MS. Since \\G' CHSH \\oo = 1, we have then that VP G Q, 
(P, G' chsh ) < 2 and, hence, (P, Gchsh) < 2yf2. 

Since the CHSH Tsirelson bound is achievable by a 
quantum behaviour arising from certain measurements 
on a maximally entangled two-qubit state and G' CHSH 
is orthogonal, this also shows that for this behaviour 
| |P|1 1 = 2, achieving the bound of Theorem 1. It is a 
natural question to ask which other behaviours in Q\C 
can attain it. We computed ||P||i for the quantum be¬ 
haviours yielding the largest known value for several two- 
outcome Bell inequalities given in [30] but, in general, 
the bound of Theorem 1 is not attained. Interestingly, 
when this occurs, the behaviour arises from a maximally 
entangled state of qubits. This seems to extend for sce¬ 
narios with more outcomes. In particular, in (2233) the 
quantum behaviour maximally violating the CGLMP in¬ 
equality [31] was given in [13, 32]. However, for it we find 
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that ||P||i ~ 1.98, while the maximal value ||P||i = 2 
is attained for the behaviour that yields the maximal 
CGLMP value when restricted to a maximally entangled 
two-qutrit state [31, 32]. This leads to consider whether 
for every bipartite maximally entangled state of dimen¬ 
sion d, 


d-l 




3 =0 


(9) 


there exists measurements such that the corresponding 
behaviour attains the bound of Theorem 1. In the fol¬ 
lowing we show that this is indeed the case 'id. We will 
use the construction of [31] that provides measurements 
on a (22 dd) scenario (i. e. x ,y = 1, 2 and a 7 b = 1,..., d) 
such that acting on |$J) lead to the behaviour 


P d (ab\xy) = (2d 3 sin 2 [n(a - b + a(x) +/3(y))/d}) 

( 10 ) 

where a(l) = 0, a(2) = 1/2, /3( 1) = 1/4 and /3(2) = 
— 1/4. This construction is enough to show that for a 
maximally entangled state of any dimension, there ex¬ 
ist a behaviour for which the bound of Theorem 1 can 
be saturated. However, it should be stressed that other 
behaviours arising from maximally entangled states can 
have this property as well. 

Theorem 6. For any maximally entangled state | < f > ( /}, 
the corresponding (22 dd) behaviour given in Eq. (10) at¬ 
tains the bound of Theorem 1: ||Pd||i = 2 id. 

Proof. In order to obtain ||Pd||i we will compute its 2d 
singular values, which we will denote by {crd~, q for 

reasons that will become clear later. To that aim we will 
first show that P d is normal (i. e. P\P d = P d P d ) implying 
that the singular values correspond to the absolute values 
of the eigenvalues. By fixing the inputs x, y 7 the matrix 
of the behaviour (10) can be partitioned as 


Pd = 




( 11 ) 


where 


where it should be understood here that d + 1 = 1. 
We will use the following properties of d x d circu- 
lant matrices [33]: they all have the same eigenvectors 
...,/\fd corresponding to eigenvalues 

d 

Xj(Pxy) = '^2P(lb\xy)u} b j ~ 1 , j = 0, l,...,d- 1, (15) 

b= 1 


where uij = exp(2nij/d) are the d-th roots of unity. This 
particularly implies that all circulant matrices are normal 
and commute with each other. Hence, it is easy to check 
that P d is normal if B = C^C. To see that this is 
indeed the case, consider that Eq. (15) tells us that the 
eigenvalues of our matrices are given by 


d-l 




0J k 

3 


2d 3 — sin- 

k—0 


( 7T k _ 7T \ 

^ d Ad) 


d-l 




. .k 


2d3 ti ™ 2 df + fi) 


k=0 

d-l 


V c > = «r £ 


2W3 q ;,,2 ( TTk _ 37T1 ' 

Za fc=0 Sm (~ 4 d) 


(16) 


These summation formulas are computed in the Ap¬ 
pendix, obtaining 


A j (A) = - l -^-(j + i(d-j)), 

Xj (B) = l -^(j-i(d-j)), 

^(C) = ^(j-i(d-j)). (17) 

Notice then that |Aj (B)\ = |Aj(C)| corresponding to the 
same eigenvector, which in addition to the fact that these 
matrices are normal, implies indeed that B = C^C. 
Hence, P d is normal and, therefore, its singular values 
are the absolute value of its eigenvalues. We compute 
now then the latter. Since all circulant matrices are di¬ 
agonalized by the same unitary U (i. e. X = UDxU 1 
with D diagonal for every circulant matrix A), it holds 
then that the matrix 


P(ab\ll) = P(ab\22) = (2d 3 sin 2 [7r(a — b + 1/4)/d]) , 

P{ab\\2) = (2d 3 sin 2 [7r(a — b — l/4)/d]) 1 , 

P(a6|12) = (2d 3 sin 2 [7r(a — b + 3/4)/d]) 1 . (12) 

Thus, to ease the notation we will rewrite Eq. (11) as 

P ‘ = ( C A ) ■ |13 > 

Notice that the d x d matrices A, B , C are circulant, as 
for every fixed x, y it holds that 



has the same eigenvalues as P d . Using now the Schur 
complement condition that tells us that 

det(.D) = det(DA) det(ZUt — DbD^ 1 Dc), (19) 

we have that 

(2o) 

and, hence, 


P d (ab\xy) = P d {a + 1, b + 1| xy) 


(14) 


A f(P d ) = A M) ± sJ\{B)\ 3 {C) 


( 21 ) 
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Thus, using Eqs. (17), we obtain that 


Oj \*d) I Aj K^d )| ^2 1 

( 22 ) 

which means that 

= ) = < 23 > 
Hence, we finally obtain that 
d-1 

m = + ^7( p rf)) = 2- (24) 

i=o 

□ 


This result is interesting because it shows that the 
bound of Theorem 1 is attainable in Q\C and that be¬ 
haviours arising from maximally entangled states are ex¬ 
tremal in this sense. Notwithstanding, it has some fur¬ 
ther application. If a real square matrix P has singular 
value decomposition given by 


P = a i U i V i = UE V T , (25) 

i 

then 


||P||i =tr£ = tr(PO) (26) 

with O = VU T orthogonal (as so are U and V). Thus, 
by choosing G = O t , we can always construct a Bell ex¬ 
pression such that for any given P it holds that (P, G) = 
||P||i. Remarkably, if we happen to have a quantum be¬ 
haviour such that ||P||i = \Jm,AfriB (i- e. it saturates the 
bound of Theorem 1), then the aforementioned prescrip¬ 
tion immediately yields a Bell expression which is maxi¬ 
mized in Q by P. This is because ||G||oo = 1 and, hence, 
by inequality (2), there cannot exist any other quantum 
behaviour R such that ( R , G) > y/mArriB■ Thus, this 
allows to construct Bell inequalities which are maximally 
violated in Q by different behaviours of interest. Theo¬ 
rem 6 shows that this is possible for maximally entangled 
states of any dimension [34]. For example, following this 
prescription for the behaviour (10) with d = 3 leads to a 
Bell expression G(4>[]’) which is then maximized in Q by 
|$+) and is given by the following coefficient matrix [46] 


G = 


/ 2+%/3 
6 

2-fa 

6 

2+V3 

6 

2-y/3 

6 

_ 1 
6 


2-V3 

2+%/3 

6 

_ 1 
6 

_ 1 
2+V3 

2-V3 

6 


_ 1 

2+\/3 

_ 1 

2-fa 

2-V3 

2 +V 3 

6 ^ 

6 

6 _ 

2 +V 3 

1 

2 -V 3 

6 

6 

6 

2—v/3 

1 

2+V3 

6 


6 

_ 1 

2-V3 

1 

2+\/3 

2+V3 

2 -fa 

6 

6 

6 


2—\/3 \ 
6 

_ 1 

2+V3 

6 

2— y^3 

2+y/3 

6 

4 / 


(27) 

Notice that this Bell inequality separates C from Q as it is 
straightforward to find that the maximal value of G(4>[j") 
under C is (3x/3 + 5)/6 ~ 1.70. A different example of a 


Bell inequality maximally violated by ]${") can be found 
in [35]. 

Thus, Theorem 6 also shows that that for a maximally 
entangled state of any dimension there always exists a 
Bell inequality that is maximally violated in Q by it and 
how to construct it [36]. In this sense, one can then 
devise a DIQIP protocol for which maximally entangled 
states are optimal within Q. This might also be of rel¬ 
evance in the context of self-testing [37] if it turned out 
that the behaviour (10) is the only one maximizing these 
Bell expressions in Q. Self-testing arises when a certain 
behaviour is the unique to attain a particular Bell value. 
This allows to check the performance of a quantum set¬ 
up without trusting any of the devices, particularly when 
it can be made robust [38, 39]. Using the techniques of 
[39] with the Bell inequality G($]|") and its generaliza¬ 
tions for other dimensions, it could be possible to check 
whether robust self-testing of maximally entangled states 
is possible in this way. 

Notice, moreover, that this observation above that al¬ 
lows to construct a Bell expression such that (P, G) = 

| |P|1 1 can be used in other contexts. For instance, the 
results of Theorems 4 and 5 imply that for every PR box 
considered there we can write down constructively a Bell 
expression, and hence a potential DIQIP protocol, whose 
performance has a quantitative gap with any quantum 
behaviour. 


V. CORRELATION SCENARIOS 

If we restrict ourselves ourselves to two-outcome sce¬ 
narios (dA = ds = 2) and taking a,b £ {—1,1}, all 
behaviours in AfS can be alternatively characterized by 
the correlators 


(A X B y ) = ]T abP(ab\xy) (28) 

ab 

and the marginal expectations 

(A x ) = Y / aP(a\x), (B y )=J2 bP (b\y ). (29) 

a b 

As mentioned before, it has been shown in [15] that for 
the particular class of correlator Bell inequalities, 

YG xy {A x By) < [[GllooVTOArna (30) 

xy 


must hold VP € Q and every real tua x ms matrix G. We 
will denote by G the mi xniB matrix with entries C xy = 
(A x B y ). It can be shown that the result above can also be 
proved using similar techniques as in Theorems 1 and 2 
[40]. For that, one just needs to consider that behaviours 
in Q 1 must fulfill [13] that a ( ttia + ms) x (tua + tub) 


real positive semidefinite matrix 


Q c\ 
C T r) 


exists with 


Qxx = 1 and R yy = 1 and proceed as in the proofs of 

































Theorems 1 and 2 to upper bound tr(GG T ). This does 
not only provide an alternative proof of inequality (30) 
but it also shows that this bound cannot give stronger 
constraints than Q 1 . Moreover, as in Theorems 1 and 2 
this leads to the following condition 

\\C\\i<^mS VP eg. (31) 


However, it turns out that this condition is strictly 
weaker than Theorem 1 (i. e. every behaviour detected 
as non-quantum by the above inequality is also non¬ 
quantum by Theorem 1) as we show in the following. 

Proposition 7. For every two-outcome behaviour P e 
MS, «/||G||i > y/rriAmB, then ||P||i > y/mAm B ■ 

Proof. Using again the mapping from P to P' as in the 
proof of Theorem 4, we obtain that 


\\P\\i = 


pf p' 
-'ll - r 12 

P' P' 
■*21 U 22 


1 


> 2 (| l E P a & lli + ll^i - P [2 - P'21 + P22 111 ) 

ab 


/i ••• 1 


\i ••• 1 

y/m A m B + He'll! 


IlCIh 


(32) 


where the inequality comes from Corollary 3 in [41]. □ 


The above result agrees with intuition since with G we 
are not fully characterizing P. This suggests to consider 
on the analogy of Theorem 2 a condition including the 
marginals (A x ) and ( B y ). Indeed, using similar argu¬ 
ments based on the first step of the NPA hierarchy one 
can show that, for every matrix G G it must 

hold for every quantum behaviour that 


y G xy ((A x By) - (A x )(By)) < \ \ G \ | yj m a m B 

xy 



(33) 


This particularly implies the following. 

Theorem 8. For every two-outcome behaviour P € Q it 
holds that 


ini, 

where C' has entries C xy = ( A x B y ) — (A x )(B y ). 

Proof. We have to show that inequality (33) is true for 
quantum behaviours. This follows from the fact that Q 1 


is equivalent [13] to the positive semi-definiteness of 

( 1 (A) (B)\ 

(A) T Q C , 

V (B) T C T R ) 


where (A) ((B)) is an tua( m B )-dimensional vector with 
entries ( A x ) (( B y )). By Schur’s complement condition 
this leads to 



with Q' xx = 1 — ( A x ) 2 and R' yy = 1 — ( B y ) 2 . Proceeding 
as in the proof of Theorems 1 and 2 we can upper bound 
tr(G'G T ) to obtain the desired result. □ 


VI. CONCLUSIONS 


We have shown that the first step of the NPA hierarchy 
allows to obtain simple analytical conditions constraining 
the set of quantum behaviours in general bipartite Bell 
scenarios, whose strength and non-triviality have been 
illustrated. Since not all problems in quantum nonlocal¬ 
ity and DIQIP can be addressed numerically, we expect 
these conditions to be of utility, filling the hitherto lack 
of such general tools. In fact, we have applied these con¬ 
ditions here to obtain a variety of results. In Sec. Ill we 
have shown that our bounds are tight and, in general, 
non-trivial, and we have used them to prove that there 
exists a finite gap (whose size we have estimated) be¬ 
tween the quantum set and PR boxes in several general 
scenarios answering a question raised in [14]. In Sec. IV 
we have provided a systematic construction of quantum 
bounds for arbitrary Bell inequalities and we have shown 
that for a maximally entangled state of any dimension 
one can obtain a behaviour that attains the bound of 
Theorem 1. Interestingly, this can be translated into a 
Bell inequality whose Tsirelson bound is reached by such 
a state. This could be applied for robust self-testing of 
maximally entangled states using the techniques devel¬ 
oped in [39]. Finally, in Sec. V we studied the particular 
case of correlation scenarios and established some links 
with the results of [15]. Several other ideas will be fur¬ 
ther investigated in the future. Given a Bell inequality, 
it would be interesting to find a procedure to find the 
best form of G in (2) to obtain its quantum upper bound 
and when it can be optimal. It is also worth further re¬ 
search to characterize which behaviours in Q\C attain 
the bound in Theorem 1: do they only arise from max¬ 
imally entangled states? It would be also desirable to 
extend this approach to the multipartite setting. Last, it 
is worth studying whether stronger analytical conditions 
as those derived here can be obtained by considering fur¬ 
ther steps of the NPA hierarchy. 
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and, therefore, 

sin(7r(fc + 9)/d) = _ 2 e l _ x2ng e~ a ^ k+e ^ n+1 ^ d . 

n =0 

(A4) 

Thus, we find that 


Appendix A: Summation formulas 

Here we prove the summation formulas used in the 
proof of Theorem 6: 


1 u 


— y_-= 

2 d 3 4 5 6 7 8 9 St, 2 (Eh 


■ jl12. 

te 2 d 


2d3 to sin2 (ir ~ is) d2 

i d— 1 k ■ —ill. 

1 uij te 2d 


E 


2d3 to sin2 (ir + is) d2 


-Y- 

9/73 ZZ 0 : 


k =0 
i -1 


UJ- 


■ ? 37r _J 
ie 


sm 


fc=0 

which follow from 

d-l 


/ 7rfe 37 t\ ^2 

V d 4d) 


O' +i(d-j )), 

-( j~i{d-j )), 
{j-i{d-j)), (Al) 


s »=S- 


^4 


fe=0 sin (7r(fc + 0)/d) 

aj -i2-Kj6/ d 

+ etz. 

(A2) 


^ g—72770^2 


To verify Eq. (A2), first notice that the geometric sum 
yields 


d -1 

E< 

n—0 


_ i27T0^77r(fc+0)/£7 


727r(fc+0)n/t7 _ ^ _Z_ / a 

2/' sin(7r(fc + Q)/d) ' 1 ' 


5 e = - 


^ g— 72770^2 


£7-1 

x £ e 

m,n=0 


£7-1 


—727T0(ra+n+l)/£7 V —i2n(j—m—n—l)k/d 


! E 

/c—0 


(A5) 


Notice that the inner sum is equal to zero, unless j — 
m — n — 1 = zd for any integer 2 for which the sum is 
equal to d. Given the values the indices take, this can 
only happen for z = 0 or z = -1, hence obtaining 


S e = - 


4 d 


_ g— 727t0^2 

£7-1 

x ]T e~ i2ne ^ m+n+1 ^ d (S m ^ n j-i + <5 m +nd-i+d) 

m,n—0 

4 d ( . , «: - i2-K(j+d)8/d ^ 

(A6) 


(1 - e -i27rfl)2 


( je -^/s + (d _ j)e - 


that leads to the desired result. 
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